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ABSTRACT 


The  impedance  characteristics  of  a  capacitive  ionos¬ 
pheric  rocket  probe  have  been  determined  as  a  function  of 
the  electron. density,  collision  frequency,  and  temperature* 
The  physical  model  used  in  developing  the  theory  incorporated 
both  effects  of  the  ion  sheath  which  forms  about  the  rocket 
body  and  the  anisotropic  nature  of  the  surrounding  magneto- 
ionic  medium.  Results  are  obtained  and  numerically  evaluated 
to  show  that  the  sheath  makes  a  very  significant  contribution 
to  the  Impedance  characteristics* 
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INTRODUCTION 


The  general  oharaeteristlee  of  the  lonoepherle 
electron  density  profile  below  120  km  have  been  determined 
by  several  methods  using  instruments  carried  by  rockets* 
These  methods  are  of  great  interest  since  ground  based 
radio  techniques  have  not  lead  to  reliable  determinations 
of  the  electron  density  profile  below  120  km*  In  1950 
some  data  were  obtained,  between  30  and  80  km,  by  the  use 
of  rocket  borne  Oerdien  capacitors  at  White  Sands,  New 
Mexico,  giving  one  of  the  first  estimates  of  the  conduc¬ 
tivity  of  this  region.  Bourdeau,  Whipple,  and  Clark  (1959) 
analysed  the  Oerdien  capacitor  data  based  on  the  assumptions 
that  cosmic  radiation  is  the  principle  souree  of  ionisation 
at  lower  altitudes,  and  that  the  contribution  of  free 
electrons  to  the  negative  conductivity  is  negligible  below 
80  km*  The  results  appeared  to  indicate  that  the  effect 
of  the  shock  front  surrounding  the  rocket  on  the  conduc¬ 
tivity  was  not  significant* 

Since  about  1957*  a  great  deal  of  work  has  been 
done  towards  developing  a  reliable  rocket  probe  applicable 
to  the  lower  portions  of  the  ionosphere.  Some  of  the  more 
popular  methods  are  radio  frequency  impedance  probes  for 
measuring  the  electron  density,  Langmuir  probes  for  electron 
temperature,  electron  density,  or  ion  density,  Oerdien 
capacitors  as  mentioned  above,  and  combinations  of  these 
and  other  methods  employed  simultaneously  to  utilise  the 
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beat  features  of  each. 

Tha  RF  proba  taohnlqua  for  the  determination 
of  ionoaphorio  parameters  uaaa  tha  inpadanoa  character** 
iatloa  of  a  proba  immersed  in  a  plasma.  In  flights  of 
RF  probea  discussed  by  Jaokaon  and  Fickar  (195?) >  tha 
ionosphere  was  observed  to  have  an  effect  on  tha  lmpedanoe 
of  a  rookat  borne  antenna  operated  at  high  frequencies. 

This  affeot  was  later  used  as  the  basis  of  an  alaotron 
density  proba  by  Kane,  Jaokaon,  and  Whale  (1962).  They 
obtained  an  electron  density  profile  which  agrees  in  form 
with  the  profile  reduced  from  lonosonde  data  above  120  km. 

In  this  paper  an  attempt  was  made  to  account  for  the  sheath 
effeots  by  considering  tha  sheath  to  be  an  empirically 
determined  cylindrical  capacitor  which  was  assumed  to  be 
ooncentrio  with  the  antenna.  The  anisotropy  of  the  lono> 
sphere  waa  alao  considered.  However  Kane,  Jackson,  and 
Whale's  analysis  is  only  applicable  where  the  collision 
frequency  is  low  compared  to  the  collision  frequency  in 
the  $0  to  120  km  region. 

An  analysis  of  a  rocket  probe  experiment,  which 
aeems  more  suitable  for  the  £0  to  120  km  region  of  the 
ionosphere,  was  made  by  R.  F.  Mlodnosky  and  0*  K.  Qarriott 
(1962).  In  this  experiment  a  small  low  frequency  RF 
voltage  was  impressed  upon  a  capacitor  built  into  the  nose 
of  a  rocket.  When  the  rocket  passed  through  the  ionosphere, 
the  impedance  of  the  capacitor  changed  in  response  to  the 


changes  in  the  electron  density,  collision  frequency,  and 
temperature.  The  frequency  of  the  voltage  was  high  enough 
so  that  the  ions  present  were  unable  to  follow  the  voltage 
changes,  and  far  enough  below  the  plasma  frequenoy  so  that 
the  eleotron  motion  was  in  phase  with  the  applied  voltage. 
Their  analysis  assumed  that  the  capacitive  coupling  be¬ 
tween  the  sheath  and  the  rocket  would  be  the  dominant 
factor  in  determining  the  impedance  characteristics  of 
the  probe.  The  presentation  given  by  Mlodnosky  and  Qarriott 
demonstrated  the  feasibility  of  this  type  of  rocket  probe. 

The  Langmuir  probe  technique  (Langmuir,  1931) 
is  based  upon  the  Boltzmann  equation  for  ionized  gases. 

When  an  electrode  is  placed  in  a  plasma,  such  as  the  iono¬ 
sphere,  and  the  current  to  it  is  measured  as  a  function 
of  the  potential  applied  to  the  electrode,  the  resulting 
current-voltage  curves  permit  the  evaluation  of  the 
eleotron  temperature  and  density.  Rocket  flights  using  a 
Langmuir  probe  have  been  employed  a  great  deal  since  19U9* 
Some  of  the  people  engaged  in  the  early  Langmuir  probe 
measurements  were  0.  Hok  (1951),  I.  0.  Johnson  (1950), 
and  J.  C.  Seddon  (1953)*  However  it  was  not  until  recently 
that  the  theory  of  such  probes  as  applied  to  the  ionosphere 
has  begun  to  be  developed.  A  paper  by  L.  0.  Smith  (1961) 
discusses  the  qualitative  features  of  this  theory  quite 
thoroughly.  He  points  out  that  the  method  is  limited  by 
the  sweep  rate  of  the  programmed  potential,  and  that  the 
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theory  foils  to  Albount  for  ths  earth's  magnetic  field  or 
for  the  ion  sheath  that  fores  about  the  rodket*  A  similar 
experiment,  dlSoillsed  by  Hoegy  end  Brace  (l96l),  using  a 
dumbbell  shaped  probe,  aooounted  for  the  sheath  and  the 
megnetio  field  ih  the  theory*  However,  here  too  the  results 
are  good  only  above  120  km,  and  even  in  this  region  there 
remain  some  large  discrepancies  in  the  measured  temperature 
to  be  explained*  The  Langmuir  probe  theory,  while  well 
developed  for  laboratory  measurements,  is  still  far  from 
being  a  well  developed  theory  when  applied  to  the  iono¬ 
sphere.  In  particular,  the  Langmuir  probe  analysis  does 
not  apply  to  the  lower  altitudes. 

2.  THE  GENERAL  STATEMENT  OF  THE  PROBLEM 

A  rocket  probe  flight  similar  to  the  one  desorlbed 
by  Mlodnosky  and  Oarrlott  is  planned  for  the  near  future  by 
the  Ionosphere  Research  Laboratory  of  the  Pennsylvania 
State  University,  The  main  purpose  of  this  experiment  will 
be  to  measure  the  electron  density  as  a  function  of  altitude 
between  50  and  120  km.  A  low  frequency  RF  probe  was  se¬ 
lected  because  of  its  ability  to  single  out  electrons  while 
rejecting  ion  and  dust  particles  on  the  basis  of  their 
respective  masses.  Therefore,  it  is  desirable  that  a 
theoretical  investigation  of  the  probe  be  conducted  to 
aid  in  the  interpretation  of  the  experimental  data. 

The  specific  problem  to  be  considered  will  be 
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the  investigation  of  the  contributions  to  the  impedanoe 
oharaeteristies  of  a  lov  frequency  EF  capacitor  probe 
immersed  in  the  ionosphere  as  a  function  of  altitude*  and 
to  find  the  relationship  between  the  measured  probe  im¬ 
pedance  and  the  electron  density* 
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3.  QUALITATIVE  FEATURES  OF  THE  THEORETICAL  MODEL 

The  capacitor  probe  la  built  into  tha  conical 
nosa  of  a  rocket  by  splitting  the  cone  into  two  sections 
separated  by  a  dieleetrio  slab.  This  split  cone  is  eon- 
neoted  by  a  wire  to  each  of  its  tvo  metallic  portions  to 
a  source  of  low  frequency  RF  voltage.  Unfortunately,  the 
conical  geometry  is  extremely  difficult  to  analyse  so  that 
a  similarly  split  cylinder  was  considered  instead;  Fig.l. 
substitution  of  a  cylinder  for  a  cone  will  be  shown  to 
have  little  effect  on  the  results.  From  the  application 
of  a  low  frequency  RF  voltage  to  the  capacitor,  we  wish 
to  determine  its  complex  impedance  in  a  slowly  changing, 
inhomogeneous,  conduoting,  and  anisotropic  medium. 

Finally,  the  electron  density  will  be  related  to  the 
measured  impedance. 

-If  an  undisturbed  body  is  placed  in  a  neutral 
plasma  composed  of  negative  electrons  ana  positive  ions, 
such  as  the  ionosphere,  then  both  ions  and  electrons  will 
contribute  to  a  current  from  the  plasma  to  the  body.  In 
the  ionosphere  the  electron  mean  thermal  velocity  is 
always  much  greater  than  the  corresponding  ion  velocity. 
This  results  in  a  net  accumulation  of  negative  charge  on 
the  body.  The  accumulating  negative  charge  reduces  the 
electron  flow  until  the  electron  and  ion  currents  are 
equal.  This  constitutes  the  equilibrium  condition.  Thus, 
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a  region  of  space  will  surround  tho  body  whioh  has  a  de¬ 
ficiency  of  olootrons  relative  to  a  neutral  plasma*  Snob 
a  region  is.  known  as  an  ion  sheath  ae  it  contains  a  nst 
positivo  charge  attributable  to  the  positive  ions.  The 
thiekness  of  the  sheath  depends  on  the  electron  and  ion 
temperatures  and  concentrations  as  well  as  on  the  velocity 
of  the  body  through  the  ionosphere. 

Since  the  definition  of  capacitance  is  based  on 
the  ratio  of  the  charge  stored  within  the  capacitor  to 
the  voltage  impressed  upon  it,  the  region  between  the  edge 
of  the  body  and  the  edge  of  the  sheath  clearly  makes  a 
contribution  to  the  capacitance  of  the  body.  This  capaci¬ 
tance  will  change  as  the  parameters  describing  the  iono¬ 
sphere  change.  Another  contribution  to  the  capacitance, 
in  the  ease  of  the  split  cylinder,  is  the  capacitance 
between  the  two  halves  separated  by  a  dielectric  insulator. 
Here  again  the  characteristics  of  the  medium  surrounding 
the  capacitor  are  Intimately  connected  with  the  value  of 
the  measured  capacitance.  The  combined  contributions  of 
these  effects  will  yield  the  impedance  and  capacitance 
as  a  function  of  the  electron  density,  temperature,  and 
collision  frequency. 


U.  THEORETICAL  CONSIDERATIONS  AND  DISCUSSION 


The  theoretical  discussion  will  be  divided  into 
several  distinct  steps.  First,  an  estimate  of  the  sheath 
thiokness  will  be  obtained  on  the  basis  of  a  simple  model. 
From  this  estimated  thickness  and  a  consideration  of  the 
velocity  distribution  of  the  ions  and  electrons  in  the 
ionosphere,  the  currents  to  the  split  cylinder  will  be 
obtained.  These  currents  will  in  turn  permit  a  compu¬ 
tation  of  the  capacitance  and  the  conductivity  associated 
with  the  sheath.  Second,  a  computation  of  the  capacitance 
and  the  conductivity  associated  with  the  oscillation  of 
the  electronic  charge  between  the  two  halves  of  the 
cylindrical  capacitor  will  be  made.  This  computation 
will  take  into  aoeount  the  effects  of  the  anlsotropio 
nature  of  the  ionosphere  based  on  a  tensor  description. 

An  os.' mate  of  the  sheath  thickness  may  be  made 
from  a  model  based  on  the  work  of  Hoegy  and  Brace  (1961) 
and  Jastrow  and  Pearse  (1957).  Assume  we  have  a  oylinder 
immersed  in  a  plasma  about  which  an  ion  sheath  has  formed. 
At  a  distance  rQ  from  the  center  of  the  cylinder  the  net 
charge  density  is  assumed  to  fall  to  zero.  In  the  actual 
case  the  net  charge  density  falls  to  zero  smoothly  but 
rapidly  as  the  distance  from  the  center  increases. 

However,  it  will  be  assumed  that  a  sharp  boundary  exists 
for  the  purposes  of  the  calculation.  The  positive  charge 


contained  within  tbs  sheath  exactly  cancel*  tbe  nagatlT# 
c barge  q  on  the  cylinder.  Then,  If  VQ  is  tbe  potential 
on  tbo  cylinder,  it  will  be  related  to  the  charge  q  by 


1  d 

V  3K 


<  »fi> 


u) 


with  tbe  boundary  conditions 


7  ■  0  when  R  £  r  , 

0 


V  ■  V  when  E  »  r  • 
o 


Than  the  expression  for  the  potential  ineldo  tbe  sheath 
le  given  by 
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B  In  R  ♦  C  , 


(2) 


where 


B 


1  fife 

InTTr  [iTio 


C  *  To  “  B  ln  ro  * 

o 

p  p 

The  charge  q  is  given  by  q  ■  •  *  h  K  e  (  r#  -  r), 

and 

h  *  the  helgbt  of  the  cylinder  in  neters, 

I  ■  the  quantity  of  charge  per  cubic  meter, 


r  -  the  sheath  radius  in  meters, 
o 

r  ■  the  cylinder  radius  in  meters, 
e  ■  1*60206  x  10-W  coulombs, 

V  ■  the  potential  in  volts, 

VQ  ■  the  potential  at  the  cylinder  surface* 
The  preceding  calculation  oan  be  used  to  estimate  part 
of  the  capacitance  and  conductance  0  (  in  amperes  per 
volt  )  of  the  cylindrical  ca'pacitor  in  the  ionosphere 
(  Hlodnosky  and  Oarriott,  1962  ).  As  long  as  the  plasma 
frequency  is  greater  than  the  oscillator  frequency,  the 
application  of  a  small  sinusoidal  voltage  will  cause  a 
current  to  flow  in  phase  with  the  applied  voltage*  Then 
we  may  write 


Q 


AI 


TV 


* 


where  AI  and  AV  are  measured  internally  at  the  generator* 
The  value  dV  ■  AT/2  is  measured  across  the  sheath  so  that 


Q 


1  dl 

7  TV  • 


The  amount  of  charge  variation  arising  from  changes  in 
the  sheath  dimension  with  the  applied  voltage  gives  rise 
to  a  component  of  current  in  phase  quadrature  with  the 
applied  voltage.  This  current  accounts  for  a  oapaoitive 
reactive  part  of  the  admittance  C,  and  is  given  by 
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«  *  i  *  . 

or  more  usefully, 

0  "  7  (t$St)U  -  7o  •  (3) 

Similarly  for  Q, 

0  "  1/2  (U^§)lv  -  vo  *  (U) 

la  ordor  to  use  equations  (3)  and  (It)  it  will  bo  necessary 
to  enter  into  a  calculation  of  the  currents  which  flow 
froa  the  ionosphere  to  the  cylinder.  The  computation  of 
the  lea  and  electron  ourrents  will  be  carried  out  sepa¬ 
rately. 

It  will  be  assuaed  that  outside  of  the  sheath 
the  presence  of  the  cylinder  is  not  felt  in  the  iono¬ 
sphere. 

lt.l  ION  CURRENTS 

Since  the  rocket  Telocity  will  normally  be  much 
greater  than  the  ion  thermal  Telocity,  the  main  contri¬ 
bution  to  the  ion  current  will  result  from  the  interception 
of  the  ions  by  the  speeding  rocket.  In  addition,  the  ion 
mass  is  so  large  that  the  ions  are  unable  to  respond  to 
the  Toltage  changes  on  the  order  of  100  ko.  This  means 
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that  tha  lea  eurrent  froa  the  loaoepbere  to  the  roeket 
la  Tirtually  lndepeadeat  of  the  rooket'e  poteatlal  with 
reepeot  to  the  ionosphere. 

.  •» 

Assume  that  all  the  loae  whleh  lateroept  the 
eheath  boundary  are  oolleotod  at  the  cylindrical  conducting 
surfaee.  If  a  aurface  element  dS  le  mowing  with  a  Telocity 
ux  along  the  z  axis,  and  the  following  quantities  are 
defined t 

n  ■  the  unit  normal  rector  to  dS, 
k  ■  Boltsmann's  constant  (1.38  x  10“2^joulee/°l), 
t  -  the  particle  Telocity  in  meters  per  second, 

T  -  the  temperature  in  degrees  Kelrin, 
m  -  the  particle  bass  in  kilograms, 

J  ■  the  current  density  in  amperes/square  meter, 
u  ■  the  rooket  Telocity  in  meters  per  second, 
then  the  number  of  particles  per  unit  rolume  with  Teloei- 
ties  between  t  and  r  ♦  dr  is 

dN  ■  £3 — J72  exp  (  -  T2a2  ,  dT  .#  (5) 

a1  a 

a*1  -  y/i  ’fc'  T7  \  and 


The  subscript  1  will  denote  the  ions  and  2  the  electrons. 


ill 


The  currant  density  J  is  then 


lx  ■/•<%*  »x  1  * 


The  limits  of  integration  are  given  by 
0-  *  u^f  *  , 

_  5  _  m  _ 

•  *  “  *  » 

-  t#2  «  . 

The  following  standard  integral*  will  he  used. 


2  2 

-  i  t! 


‘  T  dr  -  1  vr  . 

•m 

/»,.  dTx  •  f-— f— 

» J  •“*  T  dr  •  »  (  1  ♦  nt  »  t  )  , 


erf  u  ■ 


4T  * 


(6) 


(7»> 

(7h> 

(7o) 


where 


Using  tbs  definition  for  *  *•  *>»▼• 


J 


z 


.rf  •!»„) 


(8a) 


Kqnation  (8a)  is  readily  generalised.  Since  n  ie  the 
unit  noraal  to  dS ,  let  u  sake  an  angle  0  with  n  and  writs 
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He 

2\/s  aj 


a^u.n  (  J.  ♦  erf  a^u«n  ) 
-  a2  (  u»n  )21 

♦  s  1  J  . 


(8b) 


Safer ring  to  Figure  3, 


u*n  -  u(  u»i  )  ooe  *  ♦  u(  u»^  )  ein  •  , 


where  the  quantities  u,  1.,  J,  and  k  are  unit  yectore. 

Integrating  j  over  the  cylindrical  surface  S 
yields  the  total  current  through  that  surface.  To  this 
we  must  add  the  current  J  intercepted  by  the  effective 
cross  seotional  area  in  order  to  obtain  the  total  ion 
current  to  the  cylinder.  The  total  ion  current  I  is 
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I  •  v  I 
o 


[ J  |>/i  An  cos  ♦  (l  ♦  srf  Au  cob  ♦  ) 

9  9  9  1  A  2* 

-A  u  oob  *J  d#  dB  +  r^K  J  Bu  #ln  # 


2  2  2  ' 

(1  ♦  orf B  tt  Bin*)  ♦  e"B  u  8in  *  d*  ds 


n  r2K 


y/n  a^u  cob  e  (l  +  erf  a^u  cos  e) 


2  2  2 
-a^  u  cos  *1 


where 

i  *  JLjl.  , 

yfl  2a1 

A  »  a1(u.i)  , 

B  »  a^U.j)  . 

Since  the  rocket  velocity  is  very  much  greater  than  the 
mean  ion  thermal  velocity,  the  cylinder  radius  is  used 
to  estimate  the  effective  cross  section  rather  than  the 
sheath  radius.  Vhen  the  sheath  radius  is  very  much  larger 
than  the  rooket  radius,  the  forces  within  the  sheath  are 
far  too  small  to  deflect  the  ions  toward  the  cylinder. 


It  la  bov  necessary  to  compute  several  integrals. 


2a 


Ju  ooa  *  dt 

/■ 


An  eoi  t  erf (Am  eoo  t)  dt  • 
2a 

J exp(-A2u2coa2t)  dt  • 

2a 

J  u  ain  t  dt  • 


t 


1 


Bu  tint  erf(Bu  sin  t)  dt  • 

£a 

exp(-B2u2ain2t)  dt  • 


It  la  apparent  that  1^  «  1^  •  0  •  and  1^  may  ba 
evaluated  by  naing 


j  (.) - taZsl! —  f  *,in2nt  »  , 

"  V5R»  ♦  1/2)  l 


(from  Tables  of  Functions  by  Jahnke  and  Smde,  1938),  wharo 
J  (a)  ia  a  Boaaal  function.  Thia  gives 
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I  ■  I  r  h 
o 


LaV.-lV/2  [jo(i  a2.2/*) 


i  Jx  (i  A2u2/2) 

m 


♦  2,  [Jo(-1  b2'*2/2) 

♦i  Jjt-l  B2u2/2)J 


♦2n 


[.-*•'*  •  "  -2  2 


JQ(i  A  u  /2) 


+#-B2u2/2  J0(-i  B2u2/2)] 


ar2  pT  e^u  o 
-a2u2  eos2ij 


(15) 


o«  #  (l+erfUjU  oos  i)) 


Whan  the  arguments  for  tha  Basaal  functions  ara  large  and 
where  u- i_  is  taken  to  be  aero,  aquation  (1$)  becomes 

I  -  2y/n  K  r  h  2a.  u  cos  a  ♦  .  ■  v  - ■  , 

o  L  1  *!u  608  °J 


♦*  *  V 


(16) 
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Unfortunately,  under  ordinary  ''conditions  the  asymptotic 
approximation  to  the  Bosaol  funotion  nay  not  bo  made* 

U.2  BLICTROI  OURBBHT 

The  oloetrona  muat  bo  oonaidored  from  another 
viewpoint*  Haro  both  the  amall  mass  of  the  electron  and 
the  high  moan  thermal  velocity  relative  to  the  rocket  sug- 
geata  that  the  particle  trajectories  be  examined.  From 
the  time  an  electron  enters  a  region  under  the  rocket's 
influence  to  the  time  when  it  is  either  collected  or  again 
outside  of  the  region  of  influence,  the  rocket  may  be  con¬ 
sidered  stationary. 

In  examining  the  electron  particle  trajectories 
within  the  sheath,  wo  are  primarily  into rested  in  finding  . 
the  limiting  values  of  the  velocity  that  will  permit  the 
electrons  to  strike  the  cylinder  surface  3.  The  rocket 
flight  may  be  roughly  broken  into  two  parts.  First,  a 
region  corresponding  to  the  lower  altitudes  where  the 
electrons  will  undergo  many  collisions  before  striking  the 
cylinder,  and  second,  where  the  electrons  will  be  able  to 
traverse  the  distance  between  the  sheath  and  the  rocket 
without  collisions.  In  the  lower  region  the  contribution 
of  the  electron  current  to  the  impedance  will  be  over¬ 
shadowed  by  the  contribution  arising  from  the  interelee- 
trode  capacitance.  Therefore,  collisions  within  the  sheath 
will  be  neglected  for  the  computation  of  the  electron 


ourremt  during  the  entire  flight* 

Given  a  cylinder  of  radius  r  surrounded  by  a 
•heath  of  radius  r0,  let  us  assune  that  an  eleotrea  inter- 
oepts  the  sheath  edge*  at  A  in  figure  U,  and  attempt  to 
oaleulate  the  limiting  conditions  on  its  velocity  so  that 
it  will  barely  grase  the  cylinder*  Let  the  velooity  be 
deeempoeed  into  a  eomponent  tangential  to  the  sheath  v^, 
a  oomponent  in  the  radial  direction  vr,  and  a  component  in 
the  axial  direction  v  .  The  field  between  A  and  B  is 
assumed  to  be  conservative  and  repulsive  with  respect  to 
eleetrone*  As  stated  above*  no  collisions  are  assumed  to 
oocur  while  the  electrons  traverse  the  path  between  A  and  B. 

Under  these  conditions*  the  minimum  value  the 
radial  oomponent  vr  may  have  Is 

Tr  "  (-2«V*2)1/2  ’ 


Conservation  of  energy  and  angular  momentum  requires  that 


1  «2  <t; l  ♦  rl  .  r\)  -  i  .2  (tJ,  *  4  .  -  .,0  . 


*2  ro  T*  "  *2  r  T«1  • 


▼ri  Tn  *r*  the  values  of  the  radial  and  tangential 
components  of  the  velocity  at  B.  Since  vyl  -  0  at  B* 


t6  MS 


* 


ali 


(r0/r)S  -  1 


where  “  the  electron  nass  la  kg  (9*11  z  lO*^1  kg). 

The  electron  current  nay  bow  be  calculated  b y 
as earning  a  Maxwellian  distribution  of  velocities. 
Defining  a 2  ■  (n2/2kT)X^2  ,  the  eleetron  current  i  is 


*«  Vie  ^7? 


-.2(t2+t2+t2) 

▼r  e  dTrdT@dTs  » 


«  ’•  Tr 


where  the  Units  of  integration  are  found  fron 


02 


•  ,  (*eT0/*2)X/22 

<V*>2  - x 


The  sero  Unit  on  rQ  arises  fron  the  fact  that  this  eon- 
ponent  of  velocity  is  to  be  taken  as  positive  regardless 
of  the  direction  it  nay  assune. 

Define  L  ■  2*rQhHe  a2/a^^2,  and  consider  the 
following  integrals! 


dv  e 


2  2 

-a  v 
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lu3  THE  COMPLEX  IMPEDANCE  COMPONENTS 

V*  b ave  now  calculated  the  expreeaione  for  both 
the  eleotron  end  ion  currents.  It  must  be  emphasised  that 
the  above  analysis  applies  only  when  the  voltage  changes 
arm  slow  enough  so  that  the  electrons  fully  respond  in 
phase  with  the  voltage.  The  sheath  radius  and  the  po¬ 
tential  on  the  rocket  may  now  be  obtained  when  the  rooket 
is  at  equilibrium  with  the  ionosphere. 

The  condition  of  electrical  equilibrium  demands 
that  I  ■  i.  Using  this  condition  and  solving  for  the 

potential  V  yields 
0 

?o  "  ?  ln(  lesrE  Z)  •  (l8a) 

Since  the  derivative  of  equation  (2)  is  equal  to 
-q/2*r*ohv  another  expression  for  ?0  is 

~7o  •  iS;  [p2  +  ro(ln  7  • x)]  •  (18b) 

The  elimination  of  between  these  two  equations  results 
in  a  transcendental  equation  whiob  determines  the  sheath 
radius  rQ.  This  equation  is 

F  » •  flf  [i  *  3 <ta  5  -4  <1S) 
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Qrapbieal  solution  of  aquation  (19)  on  a  digital  computer 
permits  the  use  of  aquations  (3)  and  (U)  to  obtain  tbs 
expressions  for  the  capacitance  and  conductance  arising 
from  the  ion  and  electron  currents  and  the  sheath  fluctu¬ 
ations*  These  are 


C 


*e  h 
0 

ln(ro/r)  * 


1  ,kT 


I  - 


i  ♦ 


2 


¥T 

77 


ln( 


T tt:)  • 

r/r  )  J 


(20) 


(21) 


The  computation  of  equations  (20)  and  (2l) 
mill  be  oombined  with  the  results  of  the  next  eeotion  to 
give  the  total  impedance  characteristics  of  the  rocket 
oapaoitor  probe* 

lull  THE  APPLICATION  OF  MAONETOIONIC  THEORY 

The  second  part  of  the  problem  is  the  calculation 
of  the  contribution  to  the  capacitance  and  the  conductance 
from  the  electrons  oscillating  between  the  two  halves  of 
the  split  cylinder*  The  approach  will  be  based  upon  a 
development  of  the  Appleton-Hartree  equations  taking  into 
account  both  the  magnetic  field  of  the  earth  and  collisions 
with  neutral  partioles*  The  collisions  will  be  expressed 
through  a  quantity  known  as  the  collision  frequency,  that 
is,  the  effective  number  of  collisions  per  second* 
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Consider  an  electron  under  the  influence  of  an 
oscillating  electromagnetic  field  and  the  earth' a  magnetic 
field.  The  colliaiona  which  result  from  the  oacillating 
motion  eauae  a  dissipation  of  electromagnetic  energy  in 
the  fora  of  heat  and  work.  It  may  be  shown,  Budden  (1961), 
that  the  forces  from  the  oscillations  of  the  magnetic 
field  are  negligible  compared  to  the  electric  field  forces. 
The  peak  voltage  used  will  be  quite  small,  so  that  Tory 
little  heating  due  to  the  applied  field  will  result.  This 
restriction  avoids  considerations  of  power  loss  due  to 
so  called  acoustic  waves,  as  discussed  by  Kane,  Jackson, 
and  Whale  (1962). 

first  a  tensor  description  of  the  ionospheric 
medium  will  be  developed  and  applied  by  use  of  Poynting's 
theorem  to  calculate  the  power  dissipated  in  the  medium. 
Second,  an  alternate  method  will  be  employed,  which  will 
turn  out  to  be  quite  superior  to  the  use  of  Poynting's 
theorem.  However,  it  will  be  shown  that  both  descriptions 
agree  on  a  qualitative  basis.  Finally,  the  various  contri¬ 
butions  to  the  capacitance  and  conductance  will  be  combined 
in  an  equivalent  circuit  which  represents  the  actual 
rocket  probe  electrically.  The  Poynting  vector  treat¬ 
ment  will  be  seen  to  be  very  sensitive  to  the  method  of 
computing  the  electric  field  surrounding  the  rocket  and 
to  the  geometry  used.  The  second  method,  based  directly 
on  the  definition  of  capacitance,  is  easily  applicable 


to  any  body  having  cylindrical  symmetry. 

Tha  differential  equation  describing  the  motion 
of  an  eleetron  in  an  anisotropic  medium  with  collisions 

is 


qS  ♦ 


z  H) 


"4a  • 


(22) 


where 


m  *  the  mass  of  an  eleotron  in  kg 
q  -  the  charge  of  an  electron  in  coulombs 
S  -  the  electric  field  in  volts  per  meter 
UQ*  permeability  of  free  space  in  henry's/meter 
B  «  the  earth's  magnetic  field  in  ampere-turms 


per  meter 

v  ■  the  eollision  frequency  of  electrons  with 
neutrals 

«  ■  the  angular  frequency  of  the  oscillator 
!'■  the  electric  susceptibility 
P  •  the  polarisation 
c  ■  the  tensor  polarisability 
N  -  the  electron  density  per  cubic  meter 
The  following  basic  equations  will  be  used. 
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( 

9  -  e0X'  ,  (23b) 
I  ■  *0*iMt  ,  (23e) 
P  -  P0ei<#t  ,  (23d) 
P  -  Her  .  (23«) 


where  r  is  the  displacement  of  an  electron  from  some 
reference  point  designated  at  time  t  ■  0.  Equation  (22) 
in  component  form  becomes 


•Ix  •  ms  •  l*0a* gJ  ♦  H0aHys  ♦  mvx 


aB  -my  -  P0*HX*  *  *‘o#H**  + 


aB 


ms  - 


*AoeHyx  *  VV  +  ■v* 


By  using  (23d)  and  (23a),  and  the  following  definitions, 
the  above  may  be  reduced  to  matrix  form. 


i5-H  • 

(2U) 

a 

N.2 
?  • 

(2Ub) 

sue 

V 

co  , 

(24c) 
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1  .  1Z 

“I  +  I" 


Therefore,  by  using  (23a)  and  U  -  1  -  iz,  the  inverse  of 
tha  tonsor  polarizability  is 


■1 

t  m 


-  U 


i  I. 


-1  i- 


i 


-  U 


4  X. 


i  I. 


y 


-i  i. 


-  v 


l 

I 


Inverting  this  matrix  yields  the  tensor  polarisabllltgr  o 


■  X 

U(lZ-  0Z) 


V-  I 


-iUT  -II 

mXy-  XxX# 

i  x  jr 

’  I  U2-  I2 

-1VT  •  X  X 

y  *  y 

x  y  s 

I  T  itJT  -I  X 

U2«  I2 

ix  x  i  j 

8 

When  H  is  along  the  a  axis  o  becomes. 
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Using  tho  fact  that  •  *  *c  *  «  »•  haT* 


Doflno  tha  following  quantities t 


Tharofora, 


egH 


I(*.  ■■  M _ 

r  -  (i  -  iz) 


(25*) 


*T  •  S  4 


■  •<,  -  -rrfr 


■ 


V 


H  i2-  (l  -  il)6 


Thua,  we  bare 


(25b) 


(25c) 


(26) 


The  matrix  form  of  a  in  the  aquation  above  will 
be  assumed  to  ba  a  good  daaerlption  of  tba  ionosphara'a 
alaotrioal  behavior  in  tha  ragion  balow  150  km.  In  othar 
word*,  we  ara  aaauning  an  anisotropic,  inhomogeneous, 
and  dlaaipativa  medium  that  ia  slightly  parturbad  by  an 
electromagnetic  disturbance.  Fron  this  we  will  attanpt 
to  conputa  the  energy  dissipated  when  a  capacitor  is 
immersed  in  the  medium. 

Starting  with  Maxwell's  equations  and  the  follow¬ 
ing  definitions, 


curl  H 


90 

cST 


J  -  0 


i 


(27) 


ourl  B  ♦  ||  ■  0  , 

(28) 

Di  "  » 

(29a) 

Ji  *  •ij1*  * 

(29b) 

1 

Hi  *  »*ijBj  * 

(29o) 

wa  oan  obtain  an  axpraasion  for  tba  affective  value  of  *• 
Assuming  that  all  quantitias  have  tha  tine  dependence 
a  ,  we  oan  writ a  for  aquation  (27) 

owrl  B  -  |f(S  *  J J  dt)  ■  0  , 
curl  H  -  |^(D  ♦  J)  -  0  * 

In  conponent  fora,  by  using  (28a)  and  (28b),  (27) 

baooaaa 

(curl  H)t  -  [(«i3+  ^  «ij)*j]  » 

1,3  •  1,2,3. 

D,floe  Effective  « 
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^effective*!  *  *  T»  «ij)sj  "  f13*j  *  (3°* 


them 

curl  H  -  Ineffective  .  (31) 

*1$  may  nov  b®  associated  with  the  matrix 
defined  in  equation  (26).  This  association  will  enable 
the  Pointing  theorem  to  be  applied  to  the  field  quantities. 
If  we  have  a  surface  S  in  space  containing  a  volume  V, 
then  the  net  energy  V  passing  out  through  the  surface 
can  be  computed. 

W  -  -  Jd±vr(*  s  B)  dT  , 

V 


or 


V  - 


+  E*)-|t(D  ♦  D*) 


+  (H  +  H*) 


•J?(B  +  B*) 


dV  , 


where  E*  is  the  complex  conjugate  of  B.  A  time  average 
over  one  cycle  is 


* 
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U  ■ 

IT  tt 


rf  |f*i  fij 


*3  *  *1*13*3*  *  gS>13*3  '  V*3*3!|  "  * 


If  there  are  no  nagnotio  losses,  that  Is  then 


W4T  beoomea 


»«,  ■  W[(*J*13*3  '  *lf?3*3)]  w  •  (32) 

V 

Equation  (32)  gives  the  rate  of  heat  generation  in  a 
given  voluae  of  the  medium  due  to  eleetromagnetio  sources. 
The  quantities  E^  are  just  the  amplitudes  of  the  eleotrie 
field  beoause  of  the  time  averaging.  Performing  the 
operations  indicated  in  equation  (32)  we  have 

*tT  -  JtP  |b|2(*i  -  *T>  +  |*|2cos2(E,H)(Ip-  I*  ♦  IT-  I*) 

♦  i(KH  -  K*)(E*  Bj.  -  EjBgjj  dV  ,  (33) 

where  (B,H)  represents  the  angle  between  E  and  H,  and 


B  - 


*l'*  *  |*2 1 2  *  |*,|2  , 


r*  v  ..  2iXZ(Z2  +  I2+  1) 

T  "T  (?I  I*  -  !>*♦  UZ2  * 


,#  ,  .  21ZZ 

*P  "  *p  . ""  5  * 

r  1  ♦  l£ 
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Thus,  ¥ftT  is  s  real  quantity  and  it  is  greater  than  sero* 
As suns  that  the  inhomoganisty  of  the  nediua  is  suoh  that 
Kj,  Ip,  and  Kg  vary  slowly  over  the  region  strongly 
influenoed  by  the  capacitor*  In  the  ionosphere,  with 
a  capacitor  whose  largest  linear  dimension  is  on  the 
order  of  one  neter,  the  assumption  is  excellent* 


■  J2  «t  *  *?>  I  1*1*"  *  IT  <*P-  -  **> 

r  r  (3W 

•  /  |*|2<joa2(B,H)  dV  ♦  (!„-  I*)  /  UjEj  -  ijlj)  dV 


The  integrals  that  are  contained  in  equation  (33) 
reflect  the  influence  of  the  geometry.  Unfortunately, 
for  many  practical  geometries  closed  form  integration 
is  impossible*  This  is  certainly  true  for  the  conical 


rocket  capacitor* 


U.5  TWO  DIMENSIONAL  APPROXIMATION  TO  THE  ELECTRIC  FIELD 


The  amount  of  work  necessary  to  obtain  an 
accurate  numerical  integration  of  Laplace's  equation  for 
the  split  cylinder  is  too  much  in  three  dimensions* 

A  reasonable  approximation  to  the  fields  may  be  had  by 
treating  the  corresponding  two  dimensional  problem.  It 
is  well  known  that  only  a  few  two  dimensional  electro* 
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stable  probleas  say  be  rotated  to  yield  a  three  dimensional 
selatlem  to  Laplaoe'e  equation,  but  even  eo  the  quail tatire 
features  of  the  field  will  be  retained. 

The  fringing  field  produced  by  the  eylindrleal 
oapaeltor  it  of  primary  interest,  because  it  is  only 
this  field  which  contributes  to  the  heat  dissipation 
in  the  surrounding  aediua.  The  field  surrounding  the 
cylinder,  in  the  quasistatic  case,  nay  be  found  fron 


V2t  “  0  . 


(35) 


Kquation  (35)  nay  be  solved  using  the  general 
Sohears  transformation.  This  transforaation  is  governed 
by  the  differential  equation. 


(b4-  *)/* 


n  ■  one  less  than  the  nuaber  of  vertices  in 
the  polygon  foraing  the  boundary. 

C  ■  a  constant,  possibly  complex,  which  specifies 
the  orientation  of  the  polygon  with  respect 
to  a  set  of  axes  in  the  s^  plane. 

•i  -  the  looation  of  the  i  th  vertex  in  the 
>2  plane. 


b^  ■  the  interior  vertex  angle  measured  in  a 
counterclockwise  sense. 

The  boundaries  in  the  s  plane  are  first  trsns» 
formed  into  the  u  axis  of  the  t  plane,  and  then  to  the 
«  plane  where  the  stream  funotion  s  and  the  potential 
funotion  V  fora  a  rectangular  coordinate  system.  The 
Sohwars  transformation  for  this  case  is 


ds 

<TT 


C 


vr~=T 
— —  • 


Integrating,  we  obtain 

s  -  C  (VT=T  -  tmn^Vr^T)  ♦  O'  ,  (*6) 


where  C'  may  be  set  to  sero. 

Assume  that  we  are  in  the  region  where  t^O, 
then  t  -  1«  0,  Set  Vt  "-"1  ■  im,  where  m  is  positive, 
then 


tan_1vt  -  1  -  |  m 


n 

7 


m  -  1 

m  ♦  1 


When  z  «  iA,  ra  ■  1  and 


D  C 


THE  Z,t,  AND  W  TRANSFORMATION  PLANES 

FIGURE  5 
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Then  equation  (36)  becomes 


.  -  1| i  (VTTT  -  tan’VFTT)  .  (37) 


Equation  (37)  maps  tba  uppor  half  of  ths  a  plana  into 
the  upper  half  of  the  t  plane,  and  places  the  point  B 
at  t  -  1. 

In  order  to  map  the  t  plane  into  the  v  plane, 
the  transformation  is  easily  seen  to  be 


t 


.-iWV( 


(38) 


To  verify  the  transformation,  lot  t  »  |m|«  Than 


▼  In  |m|  -  -iitV  ♦  ns  , 


s  -  -2  ln|m|  , 

n 


v  -  o  . 


When  t  ■  -  |m|  ,  then 


V  ln(-|m|)  ■  -iwV  -  ns. 


V  -  -vo  , 


V 

•  •  ~  ln|»|  • 

Therefore,  the  boundaries  are  mapped  as  required.  The 
first  transformation  may  be  similarly  verified.  The 
potential  V  ranges  from  -Vq  to  0,  and  the  stream  function 
s  from  -  •  to  ♦  ®, 

For  the  transformation  to  be  useful,  ve  oombine 
the  two  equations  (37)  and  (38)  and  write  the  total 
transformation  in  real  and  imaginary  parts. 

a  -  x  ♦  iy  -  (|  ln  |  -  bl)  ♦  i£±  (a1  ♦  £  -  |),  (39) 

where  the  quantities  r,  1,  b^,  a^,  g,  and  d  are  related 
to  V  and  s  by  the  following  set  of  definitions: 
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-  [<e  -  l)2*  f^Uostltan-1^^))  , 

(UOa) 

1 

■  [(e  -  1)2+  f^Uin^tan-1^-^))  , 

(1:0b) 

exp(i4  )  cos(*  |  )  , 

(UOc) 

vo  o 

-exp(n  w  )  sin(n  1  ), 

Uod) 

¥  o  v  0 

Ai'.  b>!  ♦  .5  , 

(UOe) 
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result  is 
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.  1J2L 


There* ore , 


¥o  0 

TE  0*  ♦  Fz  * 
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V.  r 

V  W  ^T7^. 
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where 
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-*  > ( «••  *  r  •••  f 


-  *U  a 


y-  •*»  yjj. 
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Equations  (Ul)  and  (U2)  date mine  the  quantity  |s|  2  in 
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equation  (33) « 

The  angle  between  the  eleetrio  and  nagnetie  fielda 
■net  now  be  oaloulated  at  each  point  in  spaoe.  A  co¬ 
ordinate  system  with  one  axis  tangent  to  the  earth  at  the 
point. of  the  rooket  launohing  and  another  perpendicular  to 
the  earth  makes  a  convenient  reference  system.  Measuring 
all  angles  from  the  horisontal  axis  yields 


(B,H) 


+  L  -  tan 


(  B  /  Bx)  -H'. 


H*  ■  the  angle  that  the  magnetic  field  makes 
with  the  horisontal  in  radians. 

L  -  the  launch  angle  in  radians 
It  is  now  possible  to  perform  the  integration  in 
equation  (33)  numerically.  However,  one  further  simpli¬ 
fication  is  in  order.  Suppose  there  exists  a  conformal 
mapping,  such  as  a  Schwars  transformation,  between  the  a 
plane  and  the  w  plane. 

Consider  the  integral 


J  J  cos2  (E,H)  dV  ds. 


dx  + 


3V 

5y  dy» 
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vhara  •  •  x  ♦  ijr  Md  *  •  f  *  li« 

Squaring  and  adding, 

<a»)8  ♦  (a.)8  •  [<!5>2  ♦  <|I)8]  [  (a*)8  ♦  (ay)8]  , 


vta art  the  Oaaohj  Kainann  aquations. 


9T  .  da 
151  3? 

and 

3f  .  da 
Jj  ’51  * 


Since  the  mapping  la  conformal,  (E,H)  transforms  unchanged. 

Equation  (1*3)  redness  the  difficulty  of  per¬ 
forming  the  integration*  The  limits  of  integration  are 
nov  more  dlffloult  to  determine  in  the  untransforned  space. 
The  object  is  to  select  a  region  of  space  sufficiently 
large  so  that  almost  all  of  the  power  is  dissipated  within 
the  ohoaen  volume.  If  the  limits  are  selected  in  the 
transformed  space,  and  then  carried  over  to  the  untrans¬ 
forned  space,  an  appropriate  boundary  may  be  determined. 

The  most  convenient  boundary  lies  along  one  of  the  lines 
of  constant  stream  function. 

The  integration  and  computation  was  performed 
on  an  I.B.M.  707k  eemputer*  The  results  «f  this  essi- 
putation  are  displayed  in  the  next  section.  An  interesting 
feature  of  the  computation  was  that  the  last  term  of 
equation  (33)  turned  out  to  be  quite  small  oompared  with 
the  rest  of  the  terms.  This  term  can  be  rearranged  to 
resemble  a  polarisation  term  in  the  magnetolonlc  theory 
by  some  algebraic  manipulation.  Essentially,  it  is  a 
measure  of  the  anisotropy  of  the  polarisabillty  of  the 
medium.  Thus,  for  the  region  below  1$0  km  this  anisotropy 
is  a  negligbly  small  offset.  This  result  is  just  what 
is  expected  on  the  basis  of  qualitative  argument. 

When  the  magnetic  field  H  and  the  electric  field 
E  are  specified,  equation  (33)  can  be  written 
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-  -I-  *(>,»•«)  • 

The  total  power  dissipated  by  tha  flelda  from  tha 
cylindrical  eapacltor  la  approximately  2*WaT  ■  Q,  or 

y2  V2 

Q  -  2*— f(N,v,»)  "  ,  (in  joules/seeond) 

where  Vq  la  tha  peak  voltage  and  R  the  effeetive 
resistance,  Then  it  follows  that 

*  *  (1*  ,1")  •  “*> 

Iquatlon  (Mi)  datirminaa  tha  affaotive  raaistanee  of  tha 
medium  in  whloh  tha  oapaoitor  la  immersed.  Mo  Information 
results  from  this  method., about  the  oapaoitiTe  portion  of 
the  impedance. 

It  la  obvious  that  the  procedure  just  outlined 
la  far  from  an  ideal  treatment  of  the  problem.  The  greatest 
difficulty  being  that  the  answers  obtained  are  extremely 
sensitive  to  the  method  used  in  evaluating  the  effects 
of  tha  eleotric  field  near  the  cylinder.  The  two  di¬ 
mensional  field  solution  yields  only  the  grosser  features 
of  the  eleotric  field  surrounding  the  capacitor,  sinoe 
the  cylindrical  rotation  of  this  field  does  not  represent 


k9 


the'  true  three  dimensional  electric  field*  The  information 
that  it  deee  eoatala,  nanely  the  qualitative  behavior  ef 
the  leacea  daa  to  heat  dissipation,  ia  useful* 
la*  6  APPLICATION  OP  THX  DKPINITION  OF  CAPACITAKCB 

A  better  approach,  for  the  eoaputatioa  of  the 
lapedaaee  oharacteriatiee  of  the  eapaeitor,  ia  to  retura 
to  the  baeie  definition  of  oapaeitanoe  in  any  aediua* 

In  a  strict  operational  sense  this  definition  is 


C 


-  S 
7 


* 


where  q  ■  the  oharge  on  the  plates  of  the  oapaeitor  la 
eouloabs,  and  f  •  the  potential  dlffereaee  between  the 
plates  in  wolts*  The  definition  ean  be  written  in  a 
wore  useful  fora  by  expressing  q  and  T  in  terns  of 
integrals  of  the  field  quantities* 


C 


I 

1 


D*n  dS 
B*dL 


(U5) 


where  n  -  the  unit  nornal  to  the  surface  of  the 
capacitor, 

dS  ■  a  surface  eleaent  on  the  capacitor  plates, 
and  1  ■  a  path  along  which  1  is  integrated  between 

the  plates* 

Iquatlon  (U5)  beeonea 


0 


D»n  dS 


However,  n  ■  ,  where  Sg  is  ths  alsetrio  field  at  ths 

surfaee  of  the  capacitor.  Than 


C 


Siei3gJ 


IQ 


dS 


(1*6) 

( 1*7 ) 


By  an  analysis  similar  to  that  for  aquation  (33),  va  hare 


0  *  VJ  l*J  *  (*P  “  *j)  coa2(l,H)j  dS  •  (48) 

Tha  definitions  of  the  quantities  Kj  and  Kp  are  the  same 
as  in  aquation  (33). 

If  the  same  capacitor  is  evaluated  in  free 
space,  then  its  capacitance  mill  be 


C 


dS  . 


(U9) 


Since  neither  Kp  nor  Kj  varies  over  the  dimensions  of 
the  rocket. 
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81  cos2(B.H)  dS 


I 


III  da 


(50) 
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If  ill*  net*  of  tha  rooket  is  slender,  that  la,  if  lta  apex 
axel*  !■  onall,  than  by  using  tha  angular  syunetry  af  tha 
racket,  aquation  ($0)  baeomaa 

0  a  ♦  hill  <$o*2(I,hS>  ,  (51) 

0  0  o 

where 

|(l*H)lla3t  | 

4«.J(i,Bj>  «  | -  ]  oos2(E,H)  d(I,H)  .  (52) 

Whan  tha  elaotron  dansity  la  vanishingly  snail, 
and  tha  oolllaion  frequency  la  vary  high,  and  Xf 
approach  «Q,  or  C/0Q  approaches  1* 

Tha  first  tarn  in  aquation  (5l)  displays  tha 
isotropic  portion  of  tha  oapaeltor's  eharaetarlstios, 
while  tha  saoond  tarn  gives  the  orientation  sensitive 
behavior.  The  quantity  ^os2(I,H)>  is  readily  da- 
temlned  fron  a  knowledge  of  the  launch  angle  and  tela- 
netered  aspect  data  for  every  point  during  tha  rookat's 
flight.  Tha  earth's  nagnetie  field  is  assuned  to  nak* 
a  constant  angle  with  aone  arbitrary  fixed  reference 
throughout  the  flight. 

It  is  oloar  that  the  evaluation  of  aquation 
(52)  can  easily  be  nade  for  nany  different  gaonatrlas 
in  closed  fora.  In  addition,  an  evaluation  of  aquation 


(So)  my  hi  Md«  00  0  digital  emiputar  with  great  accuracy 
aad  eoaparitive  ease.  This  ia  certainly  not  tba  aaaa  with 
equation  (iala)  • 

k.7  IQVXTALXIT  OI&OUIT  BIPRK8MTATI0I 

If  both  the  devclopaents  preiantad  hare  are 
to  ba  eoaeidered  valid,  than  they  auat  ahov  the  aaaa 
qualitative  behavior  aa  a  function  of  the  ionospheric 
paraneters.  A  nunerical  eoapariaon  of  the  two  davelopaenta 
appaara  in  the  next  aaetion. 

The  quantities  ealeulatad  auat  now  ba  ape- 
olfieally  identified.  This  will  be  dona  with  the  aid 
of  a  aeheaatie  oireult  representation  of  the  actual 
physical  situation. 

As  was  stated  before,  the  eapaoitanea  and  the 
conductance  calculated  in  aquations  (20)  and  (21)  relate 
to  the  fluetuatlona  of  the  ion  sheath  surrounding  the 
vehicle.  This  nay  ba  represented  by  a  oapacitor  shunted 
by  a  resistor  between  the  surfaoe  of  the  oylinder  and  the 
sheath  edge  for  each  half.  The  oapaoltanee  ealoul,ated  in 
equation  (51)  is  a  coaplex  nuabor  which  eontains  both 
resistive  and  reaotive  terns.  Writing  C/G  in  its  real 
and  lnaginary  parts,  we  have 
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The  quantities  a  and  b  ara  daflnad  in  tarna  of  the  ra al 
and  imaginary  parts  of  Kp  and  Kj  (  written  Re(Kp) ,  Re(Kp), 
I»(lp),  and  )• 


*a(Kp) 

X*(Kp) 

Re(xT) 

I«(Kt) 


1  ♦  V 


zz 


+  X(T2  ~  z2  -  1) 

0  (iz  ♦  zs  -  i)z  +  Uzz 

IZ(T2  ♦  Z2  ♦  1) 

(IZ  ♦  Zfe  -  l)Z  ♦  Uzz 


Than  for  tha  lmpadanee  11  af  thia  complex  element*  we  have 

Z.  •  — i —  ,  m  •  2*f 

A  1«»C 

Z-  •  - . f  *  i  - - -a . —a-  • 

1  «#C  (a2  +  b2)  »C  (a2*  b2) 

o  o 


Identifying  tha  real  and  imaginary  parts  as  tha  series 
resistive  R  and  reactive  X  parts  respectively,  yields 

B  B 
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-b 

1  v- . 'w111 
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Rfl  and  Zs  may  be  converted  into  their  equivalent  parallel 
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circuit  components,  &  and  Xp  respeotivel y,  by 


I*  ♦  x2 

a  • 

- KT*“ 


kP 


»! ♦  *5 
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Voting  that  the  oondnetanee  0  In  equation  (21)  la  Juat 
the  reciprocal  of  the  resistance,  we  nay  now  complete  the 
Identification  of  the  resistive  and  reactive  components* 
Labeling  the  capacitance  derived  in  equation  (20)  aa  Cg 
and  the  resistance  from  equation  (21)  as  Rg,  the  quantities 
Rp,  Xpt  Rg,  and  Cg  may  be  combined  into  a  single  capaoitor 
Opf  and  a  single  raslator  Rpf  shunting  the  oscillator, 
the  total  lnpedanoe  of  the  raaultlng  oireult  will  ba  the 
value  aotually  measured  by  aa  Impedance  bridge  oontained 
within  the  rocket*  The  reduction  of  the  physical  oireult 
to  one  which  contains  only  one  reactive  element  and  one 
resistive  element  in  parallel  with  it  is  quite  standard* 
Kach  half  of  the  rocket  capacitor  interaets 
with  the  sheath  to  produce  a  capaoitance  Cg  and  a 
resistance  Rg*  These  four  components  can  be  reduoed  to 
their  series  equivalents  Rga  and  Xga* 


2R2 

l2*  "  i  ♦  «^c|r| 
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Ohaaging  R2>  and  Z2|  Into  two  parallel  components  Bp2 
and  Xp2  by 


and 


Combining  Rp2  and  Xp2  with  &p  and  Xp  we  obtain  the 
expresalons  for  RpJ  and  XpJ*  These  are 


(5$) 

(56) 


Xquations  (55)  and  (56)  are  the  relation 


between  the  measured  Impedance  and  the  electron  density* 


5.  IUKKS1CAL  EVALUATION  OF  THE  THEORETICAL  RESULTS 


Equations  (WO,  (55),  and  (56)  vara  numerically 
evaluated  with  ttaa  aid  of  an  I.  B.  M,  computer.  Tha  data 
for  tha  alaetron  density,  temperature,  and  eolliaion 
fraquanoy  vara  obtained  from  ionospheric  models  as  inoorpo 
rated  into  CIRA  1961.  An  experimental  determination  of 
tha  roeket's  velocity  made  from  similar  rockets  during 
test  flights  yielded  tha  approximate  result 
v*  -  3.35  X  106  -  16.85  x  103  h  , 
for  altitudes  greater  than  UO  km.  The  altitude  h  is  in 
kilometers  and  the  velocity  v  in  meters  per  second.  The 
rooket  launch  angle  vas  taken  at  85°  with  respect  to  the 
horisontal,  and  the  angle  (u,s)  between  the  rooket  axis 
and  its  velocity  vector  for  every  point  of  its  trajectory 
is  approximately  given  by 

(u,s)  ■  tan”1  10.5  -  tan“1(-2.5  x  10“^  d  ♦  10.5) 
vhere  d  -  ii.2  x  101*  (  1  i  V0,88  -  U.6  x  10“3  h 
The  data  for  the  magnetic  field  strength, 

6.07  x  10“5  vebers  per  square  meter,  and  angle  vith 
respect  to  the  horisontal,  58°,  vere  used  as  approximate 
figures  for  vhat  might  be  encountered  at  the  launch  site. 
The  oscillator  frequencies  vere  taken  to  be  1.0  x  10"*  and 
5.12  x  10^  cycles  per  second  to  correspond  to  the  planned 
rocket  flight. 

The  results  of  these  computations  are  shovn  in 
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til*  accompanying  graphs.  Figur*  6  displays  th*  ionospheric 
models  used.  Figure  7  shows  the  resistive  component  of 
th*  oomplex  impedance  as  a  function  of  altitude,  end 
similarly  Figur*  8  for  the  reaotiva  component.  Th*  values 
given  in  Figures  7  and  8  neglect  the  resistive  and  reactive 
contributions  from  the  inside  of  the  rocket.  Thus,  to 
actually  interpret  experimental  results  the  data  must  be 
adjusted  to  account  for  the  internal  contributions.  Since 
th*  interior  of  the  rocket  is  sealed  from  the  effects  of 
the  ionosphere,  the  adjustment  to  the  data  should  be  a 
constant  throughout  the  flight. 

Figure  9  is  a  plot  of  the  resistive  component 
attained  in  equation  (1*1*)  and  the  resistive  component 
obtained  from  equation  (53)*  Since  these  are  supposed 
te  represent  th*  same  circuit  elmment,  at  least  quali¬ 
tatively,  they  have  been  plotted  together  for  purposes 
of  comparison. 

The  left  hand  ordinates  of  Figure  7  to  Figure  11 

are  multiplied  by  the  capacitance  CQ  measured  in  free  space, 

while  the  right  hand  ordinates  give  the  resistance  and 

reaetanoe  for  C  -  1.0  x  10”11  farads.  The  value  of 

o  o 

selected  is  representative  for  this  type  of  experiment. 

The  application  of  magnetoionic  theory  alone, 
as  developed  in  section  l*.l*,  leads  to  probe  impedance 
curves  whose  general  form  is  not  altered  over  a  wide  range 
of  oscillator  frequencies  below  2  me.  The  frequency  2  me 
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was  an  arbitrary  computational  limit.  The  general  form 
of  tbe  reactance  curves  for  the  above  frequency  range 
is  that  shown  in  Figure  8  for  f  ■  1.0  x  10^  cps..  The 
resistance  curves  follow  the  form  of  the  curves  in 
Figure  7*  The  smooth  transition  between  curves  of  differ* 
ent  frequencies  is  the  notable  result  for  the  magnetoionie 
calculation. 

The  addition  of  the  effects  of  the  ion  and 
electron  currents  flowing  to  the  rocket  from  the  iono¬ 
sphere,  developed  in  sections  1*.0  to  k,3,  alters  the 
impedance  characteristics  considerably.  Now  the  shape 
of  the  reactance  curves  shows  a  marked  frequency  de¬ 
pendence.  At  frequencies  less  than  2.5  x  10^  cps  tbe 
reaotance  undergoes  a  change  from  negative  capacitive 
values  to  positive  inductive  values  in  the  vicinity  of 
100  km.  This  change  moves  towards  higher  altitudes  as 
the  frequency  increases.  For  the  frequency  range  2.5  x  lO'* 
cps  to  2,95  x  105  cps,  a  second  change  back  to  negative 
reactance  occurs.  Above  2.95  x  10^  cps  there  are  no  sign 
changes  in  the  reactance  until  2  me  where  the  reactance 
undergoes  two  sign  changes  again.  The  computations  were 
carried  out  in  the  range  55  km  to  130  km  because  above 
130  km  there  is  a  lack  of  meaningful  collision  frequency 
and  temperature  data.  The  curves  showing  this  frequency 
dependence  are  in  Figure  10.  Four  representative  curves 
showing  the  corresponding  resistive  component  are  in 


Figure  11 


These  results  show  the  strong  effeot  of  the 
sheath  on  the  Impedance  characteristics.  In  the  lover 
portion  of  the  altitude  range  considered,  the  oontributions 
of  the  sheath  are  only  a  small  portion  of  the  total 
impedance.  However,  at  altitudes  higher  than  approximately 
90  km  the  sheath  effect  contribute  quite  strongly.  The 
second  change  in  sign  of  the  reactance  is  directly  attribut¬ 
able  to  the  sheath.  Thus,  neglecting  the  sheath  contri¬ 
bution  or  the  m&gnetoionic  contribution  would  lead  to 
serious  error  in  the  altitude  range  considered. 

The  theory  developed  will  allow  data  for  the 
measured  value  of  the  probe  impedance  to  be  reduced  te 
electron  density  as  a  function  of  altitude.  This  is 
most  easily  done  by  computing  an  entire  family  of  impedance 
characteristic  curves  and  then  obtaining  the  best  fit  to 
the  raw  data  by  inspection,  or  by  the  use  of  a  suitable 
numerical  smoothing  technique. 
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6.  COIOLOSION 

A  relationship  between  the  measured  ippedance 
of  a  lov  frequency  RF  rocket  probe,  immersed  in  the 
ionosphere,  and  the  eleotron  density  has  been  developed. 
This  relationship  will  enable  data,  telemetered  -from 
the  rooket  to  the  ground,  to  be  reduced  into  an  eleotron 
density  profile  for  the  altitude  range  nf  $$  ka  to  130  km. 

The  problem  was  separated  into  two  main  seotions. 
First,  the  formation  of  an  ion  sheath  about  the  rooket 
probe  was  investigated  so  as  to  be  able  to  determine  its 
effeots  on  the  impedance  of  the  probe.  The  main  features 
of  the  sheath  were  determined  by  oaloulating  the  ion  and 
eleotron  emrrente  from  the  ionosphere  te  the  rooket  prebe, 
and  requiring  that  they  be  equal  at  equilibrium*  In 
order  to  apply  the  results,  the  oscillator  frequency 
must  be  low  enough  so  that  is  able  to  respond  fully  to 
the  voltage  changes.  Second,  a  magnetoionic  description 
of  the  ionosphere  was  developed  to  determine  an  effective- 
dielectric  constant  to  be  applied  to  the  capacitor.  This 
description  accounted  for  the  oscillations  of  the  electrons 
between  the  two  halves  of  the  capacitor.  Since  collisions 
were  included  in  the  description,  a  resistive  and  reaotive 
contribution  was  found  to  result  from  the  electron 
oscillations. 

A  schematic  description  of  the  electrical 
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properties  of  the  ionosphere,  ee  seen  by  an  impedance 
bridge  carried  by  the  rooket,  was  developed.  A  single 
resistive  and  reaetlve  element,  eaeh  in  parallel  with 
the  oscillator,  was  derived  from  the  detailed  sohematio 
circuit.  It  is  through  these  tvo  elements  that  the 
contributions  of  the  shedih  And  the  electron  oscillations 
vers  combined.  The  final  result  constitutes  a  relation¬ 
ship  between  the  measured  impedance  and  the  electron 
density. 

A  model  ionosphere  was  described  and  used  to 
determine  a  set  of  resistance  and  reactance  curves  for 
a  theoretical  model  of  the  rocket  probe.  These  curves 
clearly  indicate  that  the  sheath  effects  and  the 
■agnetoionie  effects  were  both  important  throughout  tbs 
altitude  range  of  interest.  Above  approximately  90  km, 
the  contributions  of  the  sheath  to  the  total  reactance 
actually  dominate  to  produce  a  double  change  in  sign 
for  certain  frequencies,  and  no  change  in  sign  for  other 
frequencies.  Alone,  magnetoionic  theory  produces  only 
a. single  change  in  sign  in  the  reactance  as  a  function 
of  altitude. 

There  are  obviously  some  drawbacks  to  this  type 
of  experiment.  In  order  to  calculate  the  impodance 
characteristics,  it  was  assumed  that  a  detailed  model  of 
the  collision  frequency,  temperature,  and  ion  mass  is 
known  accurately  as  a  function  of  altitude.  In  general, 


this  ia  not  the  oase.  Therefore,  auob  data  must  ba 
obtain  ad  from  a  suitable  experiment  flown  simultaneously 
with  tba  eapaoltor  rookat  probe,  or  froa  tbaoratieal 
aodala  of  tba  ionosphere. 

6.1  SU00S3TI0HS  FOR  FURTHER  RISSiRCH 

Soma  refinements  in  the  theory  suggest  themselves 
immediately.  The  sheath  model  should  be  obtained  from  a 
solution  of  Poisson's  equation  in  which  the  charge 
distribution  is  a  derived  function  rather  than  a  specified 
one.  The  Maxwellian  velocity  distribution  should  be 
modified  to  account  for  both  the  presence  of  a  solid 
boundary,  and  for  the  possibility  of  two  different 
temperatu rat  for  lens  and  electrons.  Fortunately,  the 

teaperature  always  enters  in  an  insensitive  way,  so  that 

* 

errors  introduced  by  an  imprecise  knowledge  of  the  temper¬ 
ature  at  the  time  of  the  flight  should  be  small.  Finally, 
the  oylinder  used  is  an  extreme  idealisation  of  the  conical 
rocket.  This  might  be  improved  upon  by  employing  a  split 
prolate  spheroid  instead. 
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